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ABSTRACT 



The associated production of J/ip + 7 at the LHC is studied within the NRQCD 
framework. The signal we focus on is the production of a J/ip and an isolated photon 
produced back-to-back, with their transverse momenta balanced. It is shown that 
even for very large values of transverse momentum (j>t ~ 50 GeV) the dominant 
contribution to this process is not fragmentation. This is because of the fact that 
fragmentation-type contributions to the cross-section come from only a qq initial state, 
which is suppressed at the LHC. We identify (/(/-initiated diagrams higher-order in a s 
which do have fragmentation-type vertices. We find, however, that the contribution of 
these diagrams is negligibly small. 
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NRQCD is an effective field theory derived from the full QCD Lagrangian by neglecting 
all states of momenta larger than a cutoff of the order of the heavy quark mass, m 
[I], and accounting for this exclusion by introducing new interactions in the effective 
Lagrangian, which are local since the excluded states are relativistic. It is then possible 
to expand the quarkonium state in terms of its Fock-components as a perturbation 
series in v (where v is the relative velocity between the heavy quarks), and in this 
expansion, the QQ states appear in either colour-singlet or colour-octet configurations. 
The colour-octet QQ state is connected to the physical state by the emission of one 
or more soft gluons by transitions which are dominantly non spin-flip and spin-flip 
transitions. Selection rules for these radiative transitions then allow us to keep track 
of the quantum numbers of the octet states, so that the production of a QQ pair in 
a octet state can be calculated and its transition to a physical singlet state can be 
specified by a non-perturbative matrix element. The cross-section for the production 
of a meson H then takes on the following factorised form: 

"(H) = £ J^iOSC^Lj)) (1) 

n={a,S,L,J} " b 

where P„'s are the short-distance coefficients and O n are local 4-fermion operators, of 
naive dimension d n , describing the long-distance physics. The short- distance coeffi- 
cients are associated with the production of a QQ pair with the colour and angular 
momentum quantum numbers indexed by n. These involve momenta of the order of 
m or larger and can be calculated in a perturbation expansion in the coupling a s (m). 
The non-perturbative long-distance factor (O^) is proportional to the probability for 
a pointlike QQ pair in the state n to form a bound state H. These matrix elements 
are universal in the sense that having extracted them in a particular process, they can 
be used to make predictions for other processes involving quarkonia. 

In fact, the importance of the colour-octet components was first noted |2| in the 
case of P-wave charmonium decays, and even in the production case the importance 
of these components was first seen [|J in the production of P-state charmonia at the 
Tevatron. The surprise was that even for the production of S'-states such as the 
J/if> or x/j', where the colour singlet components give the leading contribution in v, 
the inclusion of sub-leading octet states was seen to be necessary for phenomenological 
reasons || . While the inclusion of the colour-octet components seem to be necessitated 
by the Tevatron charmonium data, the normalisation of these data cannot be predicted 
because the long-distance matrix elements are not calculable. The data allow a linear 
combination of octet matrix-elements to be fixed || 0, and much effort has been 
made recently to understand the implications of these colour-octet channels for J/ip 
production in other processes: for example, J /if) production at the LEP || ||, the 
prediction for the polarisation of the J/ip [|IIJ, production of J /if) at fixed-target pp 
and np experiments jnj, inelastic photoproduction at HERA |12|, [13], [14j], production 



of quarkonium states at the LHC [ITH] and the predictions for the large-p^ production 
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of other charmonium resonances at the Tevatron [16[|. Recently, next-to-leading order 
calculations for quarkonium production at low-p^ have also been completed ]17| and 
make it possible to make accurate predictions for these processes. 





Figure 1: Quark tree diagrams. 
In this paper we consider the large associated production of an isolated photon 



and J/ip produced back-to-back at LHC energies. This process was first studied |18 
in the context of the colour- singlet model and more recently, the contribution of the 
colour-octet channels at the Tevatron has also been studied |L9], ^(J. The tree-level 
cross-section for the J/tp + 7 process can be obtained from the corresponding cross- 
sections for the photoproduction of J/ip which have been calculated in Ref. [T^|. It 
turns out that, because of the photon in the final state, the colour- singlet contribution 
is more important for this process than for inclusive J/ip production at the Tevatron. 
In fact, this process can be a sensitive probe of quarkonium production and it is our 
aim, in the present work, to understand how it may be used to unravel aspects of 
quarkonium production dynamics at the LHC. 



Figure 2: Gluon tree diagrams. 

The contributing subprocesses to J/ip + 7 production are 

q q - 2S+1 Lj 7, 



9 9 - 2b+ 'Lj 7 , (2) 

The Fock-components that contribute to J/ip production are the colour-singlet 3 S^ 
state and the colour-octet states 3 S\ , 1 Sq^ and 3 -Po 8 | i 2- The colour-singlet 3 S± state 
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contributes at 0(1) but the colour-octet channels all contribute higher orders in v. 
This is because the 3 S'[ 8 ' connects to the J/ip by emitting two soft gluons (both non 
spin-flip transitions and resulting in an effective 0(v 4 ) suppression) while the 
connects to the physical state via a spin-flip transition (the correct power-counting for 
which yields an an effective 0(v 3 ) suppression Q). On the other hand, the 3 -P(jfi,2 states 
connect to the J/ip by a single non spin-flip transition but since they are L = 1 states 
their production cross-section is already suppressed by 0(v 2 ), making them effectively 
ofO(t; 4 ). 
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Figure 3: Bda/dpx for J/tj) + 7 production at the LHC. 

In the a s perturbation series, we will first study all contributions to 0(aa 2 s ) i.e. 
the tree-level diagrams, and then we will study the effect of a class of higher-order 
(in a s ) diagrams which are likely to be important for the process under consideration. 
The net contribution of the various subprocess depends on various factors: 

• the initial parton flux 
1 See the erratum of Rcf. Q 
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• the order in v of the final contributing non-perturbative NRQCD matrix element, 
and 

• the pt behaviour of the subprocess. 




Figure 4: Fragmentation type triangle diagrams. 

In Figs. 1 and 2, we have shown some examples of tree-level contributions to this 
process. From our experience with large-pr J/ i> production at the Tevatron, one gleans 
the fact that is expedient to classify these diagrams in terms of the number of heavy- 
quark propagators in them. To this end we introduce some notation: the diagrams 
where only one gauge boson attaches itself to the QQ pair is called a one-vertex or 
IV diagram. Such diagrams have no heavy-quark propagators and are fragmentation- 
like diagrams. Similarly there are two-vertex (2V) and three-vertex (3V) diagrams, 
which have one and two heavy quark propagators, respectively. This classification is 
useful in understanding the systematics of J /if) production at large- j>r at the Teva- 
tron. In the case of inclusive production, it turns out that the IV diagrams (which 
are fragmentation-like) are those that dominate at large-pr- It is straightforward to 
convince oneself (by using the spin and colour projectors used in quarkonium calcu- 
lations) that the only state that can be produced in a IV diagram is the 3 Sf ] state. 
Consequently, J/if> production at large pr is completely dominated by 3 sj 8 ' production 
J7|. On the other hand, at the lower end of the pt spectrum at the Tevatron, there 
is a significant contribution from the and 3 Pq 8 \ 2 states. These come from 2V 
diagrams with the contribution from the 3V diagrams being less important. In fact, 
it is this px dependence that allows for a separate determination of the ( 3 S'{ 8 ^) matrix 
element. On the other hand, it is not possible to determine the ^S^) and ( 3 -Po 8 i2) 
matrix elements separately because they have a similar pt dependence owing to the 
fact that the production of these states proceed dominantly via 2V diagrams. 

At the LHC, the situation is somewhat different. At the large energies and large 
Pt values that will be available at the LHC, the production of J /if) is expected to be 
overwhelmed by the intermediate production of a 3 s{ 8 ' state. The contribution of the 
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Figure 5: Box digram. 



and 3 Po%2 states is very small ||15|| . The complete dominance of the fragmentation- 
type contributions for J/ip production at the LHC is also because of the fact that the 
production mechanisms in LHC which dominate are gg-initiated. This is an important 
point to keep in mind when we study J/ip + 7 production, because if we consider 
the tree diagrams for this process we find that there is no IV diagram contributing 
through a gg-initiated channel. The only IV diagram that contributes is in the qq- 
initiated channel. Consequently, we expect that fragmentation-type sub-processes are 
not as important to this process as they are for the case of J/ip production. 

To check out these expectations, we have studied the production of J/ip + 7 at the 
LHC (y^s = 14 TeV). In Fig. 3, we present the results for the cross-section Bda/dpr 
as a function of px- We have assumed —2.5 < y < 2.5 in our computations. For 
the input parton distributions, we use the MRSD-' and evolved them to the scale 
Q = Mt- For the numerical values of the relevant non-perturbative matrix elements 
we use the numbers tabulated below which have been obtained [|] by fitting to the 
CDF data. We would like to point out here that the inclusion of soft-gluon radiation 



effects [^2|, |23], |14j lead to lower fitted values for the non-perturbative parameters, but 
these effects are somewhat model-dependent. For the purposes of the present analysis, 
we prefer to use the matrix elements derived by Ref. 0. 

0(^( 3 Si)) = 1.2 GeV 3 , 



0^( 3 Si)) = (6.6 ± 2.1) x 1(T 3 GeV 3 , (3) 



" ,/ + ^ 8 y 0) < = (2.2 ± 0.5) x 10- GeV 3 . 



Since only the sum of ^Sq^} and ( 3 -P(| 8 ') matrix elements can be extracted from the 
J/ip CDF data, we make prediction by considering the maximal case, i.e. saturating 
the sum by either the or 3 P(j 8 ' matrix elements. Using heavy quark spin symmetry 
the other matrix elements are related (pl^^Pj)) = (2 J + 1) (oi /lp ( 3 Po)). First let 

us consider the case in which the ( 1 5 , q 8 ') saturates the sum in Eq. 3. The contributions 
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would hence come from the 1 Sq and the 3 5j ' terms. We find that even upto a value 
of 50 GeV in p?, the dominant contribution is that which comes from the 2V diagrams 
involving the ^jj 8 ' state. This is because of the fact that this subprocess is gg initiated. 
It is only above apr of 50 GeV that the gg-initiated IV diagram starts dominating. The 
upshot of this computation is that at the LHC energies there is an interesting interplay 
between the initial parton flux and the fragmentation-type effects, and because of the 
fact that the qq flux is small at these energies, the effects of fragmentation do not show 
up till the pt values become very large. We have also checked that the results are more 
or less unchanged even if we saturate the sum in Eq. 3 with the ( 3 Po^) matrix element. 

It is also important to check whether this is only because we have restricted to tree- 
level diagrams. In principle, it is possible that IV diagrams which are higher order in 
a s and which come from a gg initial state could modify this result. In spite of being 
higher-order in a s , being IV diagrams these are possibly enhanced by powers of pt/ttl. 
In the following we identify the possible higher-order diagrams that can contribute and 
try to estimate the magnitude of these corrections. 
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Figure 6: The ratio R of the box diagram contribution to the tree level cross-section, 
for the two cases described in the text. 

Let us consider the possible higher-order diagrams coming from a gg initial state 
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that could contribute to the signal. Since the signal that we demand is a p^-balanced 
J /if) and an isolated photon final state, the only IV diagrams that we can have are 
triangle and box diagrams with quark loops. These are shown in Fig. 4. It can 
be showed that the contribution of the IV triangle diagram is zero. The triangle is 
attached to one photon and two gluons, and vanishes due to Furry's theorem. 

The box diagram (shown in Fig. 5) corresponds to the process g g — > g 7 through 
a quark loop and the outgoing gluon fragmenting into a QQ pair which finally forms 
a J/ip through an intermediate 3 S^ state. We choose the number of quark flavour 
in the loop as four. This diagram has a non-vanishing interference with the tree level 
9 9^ QQ( 3 S^) 7. We begin by evaluating the interference term. The box diagram 
has three terms (corresponding to the three different ways that the photon can be 
attached to the internal quark lines) and three terms with the loop momenta reversed. 
It turns out that the contribution of the latter three diagrams are the same as the 
former. The individual diagrams are superficially divergent and we use dimensional 
regularisation to regulate it. Feynman parametrisation has to be done symmetrically 
and leads to a lot of simplification. This gives terms proportional to (L 2 ) n where L 
is the loop momenta and n=0,l,2. Terms proportional to n — 0, 1 are finite but the 
n = 2 term is divergent. It can be shown that on combining all the three diagrams, 
the 1/e pole cancels to give a finite part and a logarithmic term. We have extensively 
used FORM and MATHEMATICA for the calculation of the interference contribution 
and use the output of these packages directly for our numerical computations. The 
results of our computation is shown in Fig. 6, where we have plotted the ratio R of the 
magnitude of the interference term to the tree level cross-section. As before, the results 
are obtained using MRSD-' densities and a rapidity cut —2.5 < y < 2.5. The two curves 
Ri and R 2 shown in Fig. 6 correspond to the two values of the tree-level cross-sections 
obtained from saturating the sum in Eq. 3 with either of the two non-perturbative 
matrix elements. We find that the contribution of the interference is tiny, and the 
factors associated with the box diagram suppress any possible enhancement expected 
from the gluon flux and fragmentation contribution. Given that the interference term 
is so small, we expect that the box amplitude square contribution to be even further 
suppressed. 

In summary, we have studied the associated production of J/ip + 7, produced 
back-to-back, at the LHC. We find that this process gives us crucial insights into 
the dynamics of quarkonium production. In particular, we find that production via 
fragmentation-like diagrams does not dominate the cross-section upto values of 50 GeV. 
This is because of the fact that there is no such contribution in the gg-initiated channel, 
but only in the gg-initiated channel, where the corresponding parton flux is small. 
Beyond the tree level, we find that there is a box diagram with a gg initial state that 
contributes to this process, but we find that the contribution of this diagram to be 
negligibly small. 
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